Quantum Phase Transitions in Dissipative Tunnel Junctions 
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The Ueda-Guinea model of a dissipative tunnel junction is investigated. This model accounts for 
final state effects associated with single-electron tunneling. A quantum phase transition emerges, 
marking a boundary between insulating (Coulomb blockade) and conducting phases. The system is 
analyzed by large- A r techniques, self-consistent harmonic approximation, and Monte Carlo methods. 
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I. INTRODUCTION 
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Over 30 years ago, Mahani and Nozieres and 
deDominicis 5 predicted the existence of power law behav- 
ior in the absorption edges for X-ray transitions in metals. 
This phenomenon is due to the influence of a suddenly 
switched-on potential, due to the (screened) core hole, 
on the electrons. There are two aspects to this physics. 
One is the orthogonality catastrophe 6 due to the change 
of electronic wavefunctions in the presence of the core 
hole. The other is formation of an excitonic resonance 
between the liberated electron and the core hole^. Re- 
cently, several authors have noted the relevance of such 
nonequilibrium effects to mesoscopic systems such as tun- 
nel junctions^ and quantum dot o 9 : 10 : 11 . 

The model we will study has been described in the 
recent work of Bascones et alH, based on the original 
work of Ueda and Guinea^. Briefly, we consider a junc- 
tion consisting of two banks, left (L) and right (R), each 
described by a noninteracting Hamiltonian of the form 
= Ek £ a( k ) c kQ C ka' where a = L,R. The tun- 
neling between these banks is described by Wtunnei = 
—t e 1 ^ J^k k' c ikR c k'L + H.c. The phase variable is a col- 
lective coordinate conjugate to the charge transfer Q 
between left and right banks: [Q,e 1 *] = e 1 ^ . Asso- 
ciated with this charge transfer is a Coulomb energy, 
= (Q — Qoffsot) 2 /2C, where C is a capacitance and 
QoSset = ve accounts for the charge environment of the 
junction. Finally, the interaction between electrons and 
the global charge is written as&iSii^ 
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The fermionic degrees of freedom are quadratic in the 
Hamiltonian H = Hh+'HR.+'Htunnei+'HQ+'Hint, and can 
be integrated om— , leaving an effective model whose only 
dynamical degree of freedom is the phase field <fi. When 
c^kk' = 0, the resultant effective action is that obtained 
by Ben- Jacob, Mottola, and Schon (BMS) 14 . This ap- 
proach implicitly assumes that the energy level spacing 
in each of the banks is small on the scale of the charg- 
ing energy E c = e 2 /2C and temperature k B T. When the 
effects of Hint are accounted for, one obtains a modified 



L L 











x{l-cos (000 
The kernel K(s — s') is given by 
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K(s) = 
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Here, L = E c /k B T is the dimensionless inverse tempera- 
ture, and a = goa/in 2 , where j m is the high temperature 
conductance of the junction 1 ^ in units of e 2 /h. The BMS 
model is recovered for e = 0. The parameter e is a sum 
over contributions from the two banks j = L, R, with 
£j = —{8j/ir) 2 due to the orthogonality catastrophe and 
ej = (25/ 7r) — (dj/ir) 2 if excitonic effects are relevant. 

When the tunnel junction is placed in series with a 
capacitor, forming a single electron bosi^, the exter- 
nal charge environment is accounted for by a topologi- 
cal term in the action J2J, AStop = —2nii/W[4>], where 
W[</>] = [4>(L) — (f>(0)] /2tt is the winding number of the 
phase field <f>. 

The model is a version of the ubiquitous dissipative 
quantum rotor—. The e = case has been studied 
by several author:siLi£ii2i22i2£. Physical quantities are 
smoothly dependent on a and there is no phase transi- 
tion at any finite oS. For e > 0, however, a quantum 
(T = 0) phase transition is present at a critical value 
a c (e), as first noted by Kosterlita^. Here, we investigate 
this phase transition using a large- N expansion, the self- 
consistent harmonic approximation, and finally Monte 
Carlo simulations. 



II. LARGE- iV THEORY 



The large- N generalization of this problem was first 
discussed in unpublished work by Renn— Consider the 
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5[n( S ),A( S )] = Jdsl\(j£) + X(s)(n 2 (s)-qN)\ 
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+ i a Jds Jds' K{s - a') |n(s) - n(s')| 2 , (4) 
o o 

where n(s) is a real TV-component vector. The field A(s) 
serves as a Lagrange multiplier which enforces the con- 
straint n 2 (s) = qN for all s; typically we take q = 1/N, 
so the n field is of unit length. When N = 2, one 
can eliminate the constraint with the parameterization 
n = (cos <f>, sin (f)), whence one recovers the action of eqn. 

m 

In the N — > oo limit, the action becomes dominated 
by the saddle point A(s) = A, a constant function. One 
then has 



~cj 2 n + X + a[K(0)-K(u n )} 
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where w n = 2irn/L is a bosonic Matsubara frequency, 
and 



n(s) 
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The saddle point theory is thus a Gaussian theory, with 
correlation functions 



{rii(s)nj(s')) = G(s)5i 



(8) 



\uol + A + a[A-(0) - A-(w„)] 



Extremizing the free energy F = — L~ x Tr^ exp(— S[n]) 
with respect to A yields the equation G(0) = q, which is 
to be solved for X(a, L). 

In the zero temperature [L — > oo) limit, Lu n — ► lu 
becomes a continuous quantity, and K(0) — K(lo) = 
C t U 1 " 6 , with 



r(2 - e) cos(±7re) 



Straightforward analysis of the integral shows that A van- 
ishes for a > a c (e), where 
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FIG. 1: Solution to the large-iV model at finite L with e = 
0.25 (top) and e = 0.75 (bottom). A versus a for L = 8, 32, 
128, 512 (thin lines) and L — oo (thick line). a c (e = 0.25) = 
0.5619, and a c (e = 0.75) = 0.1226. 



From lim e ^o A t = 1 and C e= o = ti", we find, restor- 
ing q — 1/N, that for small e the critical point occurs 
at a c (e) = N/2it 2 e. This is to be contrasted with the 
renormalization group result of Kosterlitz 24 , a|? G (e) = 
(N — l)/2n 2 e. Note that these results agree to leading 
order in 1/N. 

Further analysis reveals the critical behavior of A (a) 
in the vicinity of the critical point: 
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(9) where j/ = max(l, e 1 — 1). 



A. Kubo Formula 

When N = 2 the symmetry group 0(2) has a single 
generator, hence there is one vector potential. For the 
0(Af)-symmetric case, a local gauge is effected by n(s) — > 
7£(s) n(s), where 72.(s) € 0(JV) is given by 



TZ(s) = exp [L4 a (s)T a ] 



(14) 
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where the T a generate the Lie algebra o(N). These 
7}N(N— 1) generators are N x N Hermitian antisymmet- 
ric tensors; they may be chosen to satisfy the normaliza- 
tion Tr(T a T b ) = NS ab . (For N = 2 the sole generator is 
T = a v .) There are therefore ^N(N — 1) vector poten- 
tials A a (s), variation with respect to which defines the 
O(N) currents, 



SS iut [y(s)n(s)] 
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= 2ian k (s)Tl l JdsK(u- s)m(s) (15) 
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and the noise current-current correlation function 



(T/ Q n ( s )4 n (0)) = 
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The conductance (in units of e 2 /h) is computed according 
to the Kubo formula™, 
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We define g(iui n ) by g a b{iu n ) = Ng(iu n ) 5 ab . We have no 
expression for the analytic continuation of the conduc- 
tance to real frequencies. As a diagnostic of any phase 
transition, we will examine the quantity g A = g(iw n= x). 
We find 
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A similar expression was proposed by Bascones et alH: 

g B = 47r 2 a f£) (cos (cf>(~L) - 0(0))) . (19) 

At the quantum critical point, we can evaluate DC 
(ui — > 0) limit of g(iuj n , ct c , T = 0). We obtain 



.9c = tt(1 -e) ctn^Tre) 
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Note that while the critical coupling a c is g-dependent, 
the critical conductance g c (e) is universal and indepen- 
dent of q. The approximations g A and g B also take on 
universal values at criticality: 
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FIG. 2: Large- N results for a c (e) (solid) compared with the 
Kosterlitz RG value af G = l/27r 2 e (dashed). Inset: criti- 
cal conductance g c (e) within the large- N approximation, with 
N = 2. 



III. SELF-CONSISTENT HARMONIC 
APPROXIMATION 



The fully nonlinear theory with action 

L 



SMs)] = J*»\{^) (23) 
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L L 

+ a jds Ids' K{s- s')[\ - cos (0(s) - <j>{s'))] 



is replaced with a trial Gaussian theory described by the 
quadratic action 



L L 
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where 
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V(s), or equivalently G(u> n ), is treated variationally, so 
we extremize by setting 



(26) 



where F = —L 1 Tr^ exp(— Sb[0]). This leads to the 
following self-consistent equation: 



-J— = \ u? n + 2a [da (1 - cosM) K(s) T(s) (27) 
o 



where 

r( s ) 



exp 
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= exp 
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cos(w„s)) G(w„) . (28) 



We iterate these equations to self-consistency, 

The dimcnsionlcss conductance g A (a,L) is plotted for 
e = 0.1 and 0.2 in fig. |3J As the dimensionless inverse 
temperature L is increased from L = 16 to L = 256, 
the curves apparently cross at a critical point, cv c (e). For 
a > a c , the conductance increases as the temperature 
is lowered, indicating a conducting phase. For a < a c 
the conductance vanishes as the temperature is lowered, 
i.e. the Coulomb gap survives. For larger values of e, 
however, a spurious first order transition preempts this 
critical behavior, as shown in fig. 0] The solution to 
the SCHA is hysteretic, and discontinuous, provided L is 
large enough. 



IV. COULOMB GAS REPRESENTATION AND 
MONTE CARLO 

Including the offset charge, the Euclidean action for 
our system is S — So + Sint + Stop , where 



So[4>] = \J ds 



ds 



(29) 



L L 

Sint[4>] = -ajds Jds' K(s - s') cos [</>(s) - <£(s')i 30 ) 



Stop = -iv[<f>(L) - 0(0)] = -2mvW[4>] , (31) 

where W[4>] = [<j>{L) — </>(0)]/27r is the winding number of 
the phase field. In the phase representation, the topolog- 
ical term represents a purely imaginary contribution to 
the action. However, in the Coulomb gas representation, 
the action remains purely real. 

We have discarded a formally divergent constant from 
our action, 



AS = L JdsK(s) , 



(32) 
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FIG. 3: Results of the self-consistent harmonic approxima- 
tion for e — 0.10 and e = 0.20, with dimensionless inverse tem- 
perature L = 16, 32, 64, 128, 256. The crossing of the curves 
at a single point indicates a second order phase transition. 
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FIG. 4: Results of the self-consistent harmonic approxima- 
tion for e = 0.30 and e = 0.40, with dimensionless inverse 
temperature L — 16, 32, 64, 128, 256. For high temperatures, 
the curves seem to cross at a single point. However, a spurious 
first order transition intervenes at low T (large L). 
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which may be rendered finite through an ultraviolet reg- 
ularization of K(s), viz. 



ir ( \ ) n / 77 s 



1 - e 



-L sin(7r | s\ / L) / irb 



2-e 



so that K b (0) = &~ (1 ~ e) . 

The partition function for our problem is given by 



(33) 
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where the average is with respect to the bare action Sq. 
This allows us to transform the complex action (due to 
the topological term) in the phase representation to a 
purely real action in a dipole gas representation^. The 
coordinates {s^} are interpreted as locations of positive 
and negative charges. Each factor of (— S- m t) introduces 
another pair of such charges, i.e. a dipole. We now write 
the field <f>(s) as a sum of winding plus a periodic pieces: 



2ttWs 
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where 77 is a constant and W is the winding number. The 
prime on the sum denotes exclusion of the n — term, 
which is accounted for by 77. The bare action is then 



So — 



n 2 W 2 



L ' 4 
Thus, (|^K)| 2 ) = 2/<4, and 

C(s) = i([^( s )-^(0)] 2 ) o 

L 1 — cos(2nns / L) 
7r 2 ' n 2 

n=l 

= U(l- 



(38) 



(39) 



where the subscript indicates that the expression is to 
be periodically extended from its value on the interval 
s G [0, L], Summing over the winding number W and 
averaging over <p(s), we obtain the iV-dipole pair Boltz- 
mann weight, 
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FIG. 5: Single pair potentials Vc(s) and Vc{s) + Vd(s) for 
v = (solid), = 0.25 (dashed), and z/ = 0.5 (dot-dashed). 
Upper panels have L = 5; lower panels have L — 25. In all 
cases b = 0.2. 
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where P is the total dipole moment, 

JV 

^ = E(-i — 7)- 

and $3 (z I t) is the Jacobi theta- function 29 
We now use the ^-function identity, 
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to obtain 
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FIG. 6: Coulomb gas Monte Carlo results for conductance 
3a (a) at L = 10 (filled triangles), L = 20 (open squares), and 
L — 40 (filled circles) for e = 0.30 (left panel) and e = 0.50 
(right panel). A critical point marks the boundary between 
a high a conducting phase and a low a insulating (Coulomb 
blockade) phase. 



and where the location of each charge is restricted to the 
interval S [0, L], so that P E [-NL, NL}. 

Interpreting (|43|l for v = 0, we have that the system 
consists of a gas of N dipole pairs. Each of the 2N 
charges interacts with every other charge according to a 
one-dimensional Coulomb potential —aa'\s — s'\. In addi- 
tion, each charge is logarithmically attracted to its unique 
mate - the other member of its dipole pair - via a poten- 
tial V d (s) = -lnif(s) ~ (2- e)ln|s|. For v ^ there is 
an additional electric field of strength 2v present, a is the 
dipole fugacity; a — ► will suppress the appearance of 
dipoles. Hence, in the absence of any short-time cutoff 6, 
short-distance dipole pairs will proliferate without limit. 
A simple estimation of the dipole density, ignoring the 
Coulomb interactions, is ndi P ~ 6~ < - 1 ~ e V(l — e )- I n our 
Monte Carlo simulations we have used 6=1; the essential 
physics is rather weakly dependent on b, though. Grand 
canonical averages are to be computed in the usual way, 
i.e. (A) = Tr(gA)/ Tr g, where the trace entails a sum 
over all possible numbers N of dipole pairs and integra- 
tion over their 2N charge coordinates. 

In figure |S] we plot the Coulomb energy for a single 
dipole pair, 



V c (s) 



In 
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Hl(0) 
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(46) 



as well as the total dipole energy V c (s) + V d (s) for v = 0, 



v = 4, and v — | for L — 5 and L = 25. Note how 
the interaction becomes asymmetric (yet still properly 
periodic) when v is neither integer nor half odd integer. 
Furthermore, the Coulomb interaction is effectively can- 
celled when v 
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The XY phase correlation function, 



G(s, 
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is simply related to the dipole separation correlation 
function, 
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To see this, define the quantity 
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One then has 



1 00 r 

>K(s)G(s) = — Y. drQ N (T + s,r) 
a N=oi 



h(s) 



(50) 



What is computed in the Coulomb gas Monte Carlo 
calculation is the dipole separation correlation function 
h(s), from which the XY phase correlator is obtained via 
G(s) = h(s)/aK(s). The conductance is then obtained 
using l(T%|) . 

The renormalized charging energy, E* , is given by 
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Note that when the dipole fugacity a vanishes, there are 
no pairs at all and one obtains 



E r . 



, LH»(0) 
2 H L (Q) 



(52) 



which also follows from an analysis of the noninteracting 
Hamiltonian 7Y Q= o = E c (n + v) 2 . In this case the renor- 
malized charging energy E* interpolates between its low 
temperature value of E* — E c and its high temperature 
limit of E* = 0. 



A. Results 

Typical raw data for <? A (a, L) are shown in fig. [fJl Two 
phases are identified for every positive value of e: a small 
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FIG. 7: Critical coupling a c (e) (top) and critical conductance 
g A (e) (bottom), comparing data from Monte Carlo (points) 
and large- iV (smooth curve) calculations. The dashed line in 
the top figure is the Kosterlitz RG result a^ G = l/27r 2 e. 



a insulating phase in which the Coulomb gap persists 
and g A decreases as T — > (L — > oo), and a large a con- 
ducting phase in which g A diverges as T — > 0, indicating a 
nonlinear 7—1/ relation at T = 0. The boundary between 
these phases is marked by a critical value a c (e), at which 
there is a quantum phase transition. At a = a c , the con- 
ductance becomes temperature-independent (provided T 
is low enough that the system is in the scaling regime). 
Monte Carlo results for a c (e) and <?a( c ) are presented in 
fig. [7| along with comparisons to large- iV results. 

Our results for a c {e) differ significantly from those of 
Bascones et alH, who obtained a c through analysis of 
the renormalized charging energy, assuming E*(a) ~ 
[oi c — a) 1 ^. (Our values for a c (e) are approximately 
four times smaller throughout the range < e < i.) 
In fact, we find that our raw data for E*/E c versus a 
are in good agreement with those of Bascones et alM>. 
This agreement is noteworthy since their Monte Carlo 
was carried out in the phase representation whereas our 
is in the Coulomb gas representation. However, it is very 
difficult to reliably extract a c from the charging energy 
data, as our results shown in fig. [S]show. At a temper- 
ature corresponding to L = 50 (ref. Il2|) or L = 40 (this 
work), there is no detectable signature of the phase tran- 
sition at a = a c . Similar behavior is found in the large- N 
results of fig. ^ in the large-A^ theory, X(a) plays the 
role of an energy gap, similar to E*. As is evident from 
fig. n even at relatively low temperatures of 1/L w 1/32, 
extrapolation of a c based on the zero-temperature criti- 
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FIG. 8: Dimensionless renormalized charging energy E*/E c 
versus a at L = 10 (filled triangles), L = 20 (open squares), 
and L = 40 (filled circles) for e = 0.50 (top panel) and e — 0.80 
(bottom panel). Arrows indicate the location of the phase 
transition as obtained from crossing of conductance curves. 

cal ( i.e. power law) behavior is problematic. Indeed, for 
fixed 1/L, varying a takes the system through renormal- 
ized classical, quantum critical, and quantum disordered 
regimes™, and the single parameter L — 00 behavior of 
A(a) is insufficient to extract a<2&. 

Finally, we plot conductance versus charge offset v for 
e = 0.2 and e = 0.8 in figs. I§1 and [Till respectively. Two 
values of a on either side of a c are chosen, corresponding 
to opposite temperature dependences at v — 0. We find 
a curious double peak structure in the vicinity of v = ^ 
at higher temperatures, but which disappears as T — > 0. 
For v — i and e = 0.8 the conducting state prevails even 
at small values of a. This feature is emphasized in fig. 
1111 where g A (a,v — |) is contrasted for e = 0.8 and 
e = 0.1. For e = 0.8, a conducting state is observed down 
to values of a as small as 0.01, well below the v = 
critical value of a c ~ 0.07. Hence, it is possible for the 
junction to exhibit opposite temperature dependences in 
the troughs (y m 0) and peaks (y w |) of the conduc- 
tance as the gate voltage is varied. This result is to be 
contrasted with the behavior at e = 0.1, where the con- 
ductance at v = 0.5 is very weakly temperature depen- 
dent. Indeed, at e = and v — i, the inverse charging 
energy is known to diverge very weaklyS. While the ex- 
periments of Joyez et alM> seem to be perfectly consistent 
with the more familiar e = behavior, anomalous tem- 
perature dependence in the troughs has been observed in 
the conductance of quantum dots— , although it is hardly 
clear that the nonequilibrium effects which we consider 
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FIG. 9: Conductance g A versus charge offset v at e = 0.2 for 
a — 0.2 (top panel) and a — 0.5 (bottom panel) at dimen- 
sionless inverse temperatures L = 10 (filled triangles), L = 20 
(open squares), and L = 40 (filled circles). The critical cou- 
pling for v = is a c ~ 0.44. 



FIG. 10: Conductance g A versus charge offset v at e = 0.8 
for a — 0.025 (top panel) and a = 0.10 (bottom panel) at 
dimensionless inverse temperatures L = 10 (filled triangles), 
L = 20 (open squares), and L = 40 (filled circles). The critical 
coupling for v — is a c ~ 0.07. 



are uniquely responsible for this phenomenon. 
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V. CONCLUSIONS 

Nonequilibrium shakeup processes have the potential 
to drastically affect the physics of tunnel junction be- 
havior, by allowing for a conducting phase in which the 
Coulomb blockade is completely suppressed. The phase 
transition is made manifest within a large-iV approach, 
although it was first predicted over 25 years ago by 
Kosterlita24 in his renormalization group studies of O(N) 
classical spin chains with long-ranged (power law) inter- 
actions. 

The phase transition has two principal signals. First, 
the renormalized charging energy E* vanishes for a > a c . 
Second, while the conductance vanishes for a < a c and 
diverges for a > a c , precisely at the transition g c = g(a c ) 
is finite and universal (although e-dependent). Extract- 
ing the critical value a c from numerical data at finite 
temperature can be tricky, we have found. In particular, 
the renormalized charging energy shows little signal at a c 
even at dimensionless inverse temperatures on the order 
of L = E c /k B T ~ 40; this behavior is borne out explicitly 
in our large-iV studies. It is more reliable to obtain cv c 
from the crossing of the conductance curves g(a, L). 

An extension of the model investigated here to the case 
of granular systems has been recently considered^. 
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VII. APPENDIX: ABSENCE OF 
LONG-RANGED ORDER FOR e = 



Applying Mermin's classic Bogoliubov inequality 
arguments^, Simanck 38 has proven the absence of long- 
ranged order in the 1/n 2 XY chain. Here we (trivially) 
extend Simanek's work to the continuum, and show how 
the e = case is marginal, i.e. absence of order cannot 
be so proven for e > 0. 

We begin with an action 



S[P(8),4>(8)] = ds 



l P 2( s ) + l(d s 0) 2 -hco S( f>( S ) 



L L 

+ ajds Jds'K(s-s')(l-cos[<j>(a)-^(s')]) .(53) 
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FIG. 11: Conductance versus coupling a at charge offset 
v = ~ for e = 0.80 (top) and e = 0.1 (bottom), for L = 5 (open 
down-pointing triangles), L = 10 (filled triangles), and L = 20 
(open squares). For reference, at v = 0, a c (e = 0.8) ~ 0.07 
and Q c (e = 0.1) ~ 0.66. 



The Poisson bracket of A and B is defined by 

L 



{A,B}= fds 



( SA SB 6A SB 



\SP(s) 6</>(s) S(/)(s) SP(s) 



(54) 



and the average of a functional A[P(s), 4>(s)) is 
$VP $V(f>A[P,(t>]e- s ^ 



(A) 



JVPfV<f)e- s ^ 



(55) 



Note that when A is a functional of (f>(s) alone, the func- 
tional integral over the momentum field P(s) cancels be- 
tween numerator and denominator, i.e. we recover the 
modified BMS model. 

The Bogoliubov inequality guarantees^! 



(A* A) > 



i({c,^»r 



(56) 



for an arbitrary functional C[P,4>]. Following Simanek, 
we take 



A[P,4>] = Jds sin <j){s)e" 
o 

L 

C[P,<I>] = JdsP{s)e- 1 ^ 



(57) 



where uj n = 2im/L is a Matsubara frequency. In the 
L — > oo limit, the Bogoliubov inequality then establishes 
the following identity: 



dio m 2 

2n hm + \w 2 + 2aC e \u}\ 1 -< 



< 1 , 



(58) 



where m = (cos </>(0)) is the average magnetization den- 
sity. When e — 0, the integral diverges as — m 2 \n(hm), 
hence the inequality demands m — when h = 0, i.e. no 
long-ranged order. However, the case e = is marginal, 
and for any e > the integral is convergent when h = 0. 
While this does not rigorously prove the existence of an 
ordered phase for e > 0, it is at least consistent with our 
findings. 
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